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^ ' Abstract 
(N 

^"'j We discuss a simple procedure for computing one-loop quantum energies 

C/3 ' of any static field configuration that depends non-trivially on only a single 

C > ■ spatial coordinate. We specifically focus on domain wall-type field configu- 

rations that connect two distinct minima of the effective potential, and may 
or may not be the solutions of classical field equations. We avoid the con- 
ventional summation of zero-point energies, and instead exploit the relation 
between functional determinants and solutions of associated differential equa- 
tions. This approach allows ultraviolet divergences to be easily isolated and 
extracted using any convenient regularization scheme. Two examples are 
considered: two-dimensional cj) 4 theory, and three-dimensional scalar electro- 
dynamics with spontaneous symmetry breaking at the one-loop level. 
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I. INTRODUCTION 



a: 

The calculation of the one-loop energy of an arbitrary static background field configura- 
tion is a problem in field theory that has received much attention over the years. (For a small 
sampling of previous papers, see Refs. [[L]-[5[] and references therein.) Various complications 
arise when one deals with field configurations with non-trivial spatial dependence, as extra 
care is required to relate renormalization of the effective action for such a field configuration 
to the renormalization used in the perturbative sector of a theory. 

To illustrate this situation, let us recall the result for the kink mass in 1+1-dimensional 
<ft A scalar field theory The typical way of computing one- loop corrections involves the 
summation over zero-point energies computed in the presence of a kink, and subsequent 
subtraction of the contribution from the vacuum sector. The necessary regularization is a 
rather delicate procedure, and it has been recognized some time ago that the naive use of a 
momentum-cutoff regulator leads to an incorrect answer Of course, the fact that the 

one-loop calculation of the kink mass could be performed analytically is directly related to 
the special form of the kink solution for which the spectrum of the fluctuation operator in 
the kink background is exactly known. 

For a generic background field configuration, the situation becomes more problematic 
as there is no general procedure for solving the eigenvalue problem analytically, and thus 
one must rely on numerical methods. This means that the sum over zero-point energies 
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must be regularized properly, and this regularization procedure must leave one with a well- 
defined regulator-independent integral accessible to numerical evaluation. The recent paper 
0] gives one approach for providing an unambiguous regularization of the sum over zero-point 
energies, based on the subtraction of successive Born approximations for scattering phase 
shifts. Having rendered the zero-point sum finite by this subtraction, one may then add 
the subtracted terms back in a regularized form, and observe that the regulator dependence 
cancels appropriate counterterm contributions contained in the classical expression for the 
energy. 

However, the approach of Ref. || can become complicated as one goes to higher di- 
mensions (as higher order approximations for the phase shifts are needed), or encounter 
fluctuation operators that couple multiple fluctuating fields, as one inevitably does in gauge 
theories. It is therefore desirable to have a method that would allow a more computationally 
straightforward treatment of ultraviolet divergences. 

We will discuss an alternative method, specifically applicable to kink-type field configura- 
tions that depend non-trivially on only a single spatial coordinate, which naturally simplifies 
the treatment of ultraviolet divergences and also allows one to deal with the other problems 
that appear in the treatment of gauge theories. 

The crucial feature of the method we describe is the departure from the conventional 
summation over zero-point energies. Instead, we go back to the functional determinants that 
enter the one-loop effective action, and express them in terms of the solutions to associated 
differential equations with certain boundary conditions. We examine two specific examples 
to illustrate the main ideas. One natural choice is, of course, the classic case of 1+1- 
dimensional scalar 4 field theory. This simple case provides us with a consistency check, 
as the value of kink mass can be computed independently by an analytic calculation. In 
addition, working with two-dimensional 4 theory presents an opportunity to illustrate the 
relevant ideas in a simple context. 

We then examine three-dimensional scalar electrodynamics. This theory has a first order 
phase transition (as the scalar mass is varied) separating Coulomb and Higgs phases. But 
the existence of a first order transition, and hence the presence of a degenerate ground 
state, is only apparent when one-loop quantum corrections are included. Hence, computing 
the interface tension (or domain wall energy density) at the phase transition requires the 
evaluation of the effective action of field configurations which are not classical solutions of 
the underlying classical field theory. To our knowledge, our approach provides the first 
practical technique for computing one-loop energies of arbitrary Higgs field profiles with 
one-dimensional spatial dependence in a gauge theory. 

Three-dimensional scalar QED can be viewed as a low energy effective theory for four- 
dimensional scalar QED at finite temperature. The domain wall surface tension in the 
three-dimensional theory translates into the surface tension of an interface that separates 
hot and cold phases of the finite temperature theory. This quantity has been a subject 
of several investigations in electroweak theory and scalar electrodynamics has often 

been used as a toy model [|Tl],[n]]. The desire to have a way of correctly computing the full 
one-loop interface tension in electroweak theory was part of the motivation for this work. 

This paper is organized as follows. In the next section, we consider two-dimensional scalar 
4 theory. The results of this section are not new, but this toy model provides an useful 
playground for setting up notation, introducing the basic ideas of our method, identifying 
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possible problems, and developing appropriate tools for solving them. 

In section 3, we consider three-dimensional scalar QED. We show how to deal with the 
difficulties associated with the absence of degeneracy of the ground state at the tree level. 
We also show how the scheme used to compute one- dimensional determinants in section 2 
may be extended to deal with a fluctuation operator that couples together more than one 
field. Numerical results are presented at the end of this section, where we show that they are 
consistent with expected behavior in the limit of a small ratio of scalar and gauge couplings. 

II. TWO-DIMENSIONAL c/> 4 THEORY 

In this section we will be concerned with our toy model, double well scalar 4 field the- 
ory in two dimensions. We start by formulating the theory and introducing our notation. 
We then review the traditional way of calculating one-loop corrections to the energy of 
an arbitrary field configuration^ by summing over zero-point energies. We implement this 
prescription and compute the one-loop correction to the kink mass using dimensional regu- 
larization. This warm-up computation, which reproduces old results, will serve to illustrate 
the difficulties that one faces in generalizing this approach to arbitrary field configurations 
and more complicated theories. We will argue that the integration over one of the direc- 
tions (Euclidean time), which reduces the functional determinant to a sum over zero-point 
energies, does not come without a price. As we will see, this approach forces one to resort 
to computationally inconvenient ways of regulating ultraviolet divergences. 

We then introduce an alternative method that does not rely on the summation of zero- 
point energies. We show how one can make the issue of ultraviolet regularization simple 
by integrating over the momentum in the direction of spatial variation of the field profile, 
instead of integrating over the Euclidean frequency. We relate the resulting one-dimensional 
functional determinants to solutions of associated differential equations, and present an 
explicit formula for the mass of an arbitrary field configuration that connects the two minima 
of the effective potential. The value of the one-loop correction to the kink mass is reproduced 
as a check of this general result. 

A. One-loop energy density 

Consider scalar 4 -theory with the double well potential 

wH-^+^y, (2.i) 

which exhibits spontaneous symmetry breaking at the tree level. This feature persists 
quantum-mechanically due to the — > —(f) symmetry. 

We consider a Euclidean version of this theory formulated inside a finite n-dimensional 
box with Dirichlet boundary conditions for the fluctuating fields. Formulating a theory in a 



1 In 1+1 dimensions what would be the "energy density" of a domain wall in higher dimensions becomes 
just the energy, or mass, of the kink-like field configuration. 
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variable number of dimensions allows us to use dimensional regularization, even though the 
final results of this section will only apply to the two-dimensional theory. 

We will be dealing with field configurations that vary only along a single direction. We 
will assign the subscript _L to this direction, and the relevant coordinate will be denoted by 
x±. The coordinates along all other directions are denoted by x\\. We will refer to these 
directions as "worldvolume" coordinates in obvious analogy with D-branes in string theory. 

We will denote by L the half-length of the box in the direction parameterized by x±, and 
by Vn the volume of an (n— l)-dimensional slice x± = const. Both L and V» will be taken to 
infinity at the end of the calculation, leaving us with the theory in infinite Euclidean space. 

We further restrict our attention to field configurations h(x±) that obey the following 
boundary conditions: 



h(-L) 
h(+L) 



-v, 
+v, 



(2.2) 
(2.3) 



where ±v = ((f)) are the minima of the effective potential. The tree-level value of v is given 



by v. 



2 

tree 



6 m; 



-. We require the field profile to approach the limiting values at x± = ±oo fast 



enough so that the energy density a, which is related to the effective action T[h] as 

T[h] 



a 



lim 



(2.4) 



is finite.Q 

The effective action can be evaluated perturbatively and is equal to the classical ac- 
tion evaluated at the given field profile plus quantum corrections. Retaining only one-loop 
corrections, one has the following expression for a: 



a 



dx± [C(h(x±)) 



C(v)]+ lim 

Vn,L-*c 



1 , , A 



(2.5) 



where A = — d\ — &» + V" is the fluctuation operator in the non-zero background field sector, 
while A^) is the corresponding fluctuation operator in the vacuum sector. 

We can split a into three parts that correspond, respectively, to the piece of the classical 
energy density that contains only renormalized parameters m and A, the one-loop quantum 
corrections, and the counterterm contribution to the classical energy density, 



a 



0"cl 



a, 



qu 



Or 



(2.6) 



The classical piece a c \, and the quantum correction combined with the counterterm contri- 
bution, cr qu + cr c t ., will each be UV- finite. 
The terms that enter (|2.6D are given by 



dx 



'-(d ± h(x ± )f 



nr , 2 A , , 4 3m 



2 -h(x±y + ^h(x ± y + 2A 



(2.7) 



2 We will always normalize the effective action so that it vanishes for vacuum field configurations: 
T[h = ±v] = 0. 
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(J, 



qu 



lim 



2VTi 



-tr In 



-<9| - <9 2 + V" 

-a? - a? + v(°y 



and 



Cat. 



m5 — m 



dx 



h{x x f - 



Qm 



(2.8) 



(2.9) 



In Eq. ( j2.8| ), y(°)" = 2m 2 denotes the tree-level curvature of the effective potential at 
the stationary points. 

The fluctuation operators are rcy-independent and can therefore be partially diagonalized 
by a Fourier transform over the xu coordinate. Taking the Vji —>■ oo limit recasts <r qu as 



cr, 



qu 



lim 

L— >oo 



d 



n—l ) 



(27T)™- 1 



In det 



kf\ - d\ + V" 

*f - 9i + v(°y 



(2.10) 



where detj_ is a functional determinant of a linear operator acting on the Hilbert space TIl 
of functions of x± that vanish at ±L. 



B. Summing over zero point energies 



The conventional approach to computing the one-loop energy relies on making ( |2.10| ) 
into a sum of zero-point energies. To accomplish this, one chooses one of the worldvolume 
directions to play the role of Euclidean time, and integrates over the corresponding com- 



ponent of fen 
-d 2 , + V", 



This allows one to rewrite <r qu in terms of the eigenvalues e~ of the operator 



' qu 



d 



n-2j 



(2vr) n - 



E 



(0) 



V2 



(2.11) 



This expression requires the presence of some regulator which will render finite both the 
sum and the integration over the remaining worldvolume directions. In the two-dimensional 
case, if we don't use dimensional regularization, then there are no integrations left to do and 
the expression for a qu becomes just a suitably regularized sum over zero-point energies, 



a, 



qu 



.(0) 



(2.12) 



It is expression (|2.12| ) that has conventionally been used for computing the one-loop 
energies in (l+l)-dimensional scalar field theories |],|]]. Below we will re-evaluate the one- 
loop correction ( |2.11| ) to the mass of a kink in the scalar (l+l)-dimensional 4 theory using 
dimensional regularization. The answer has been known for quite a long time Jl|, but it is 
interesting to see how dimensional regularization can naturally avoid various regularization- 
related subtleties which appear in more traditional treatments 0. 

The primary purpose of this computation, however, is to illustrate that the direct use of 
dimensional regularization applied to ( 2.11|) (without prior subtraction of Born approxima- 
tions as in H) for an arbitrary field configuration is not desirable. The computation that 
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we are about to present is feasible only because the fluctuation operator for this particular 
field profile is nice enough that its spectrum can be obtained analytically. 

Let us again go back to the theory formulated in n dimensions, keeping in mind that 
physical results will be obtained via the limit n — > 2. The quartic coupling is only multi- 
plicatively renormalized to one-loop order: 

A = /i 2 -"A, (2.13) 

with A bearing the dimension of mass squared; its beta function vanishes in the physical 
limit n — > 2. 

The relation between the /^-dependent renormalized parameter m(/i) and the bare pa- 
rameter m can be fixed by a calculation of the one-loop correction to the self-energy in 
the perturbative sector of a theory (which is defined with respect to one of the degenerate 
ground states). This correction has a pole at n = 2. To render the self-energy finite, the 
counterterm must also have a pole with the same residue. One has, of course, the free- 
dom to add some finite quantity to the counterterm. We exploit this freedom by choosing 
MS'-subtraction, which in the two-dimensional theory is defined as 

An / p 7 \ n /2-i 1 
n-2 A I e \ 1 



m ^ = m (u) 2 -u n - 2 -^ — — , (2.14) 

o v Sn\4nJ n/2-1' V 1 

where 7 is Euler's constant. 

It is easy to determine the one-loop beta function of m(fi) 2 : 

<V"M 2 = S + °< A2 >- (2 - 15) 

To calculate the kink mass we set the background field configuration h(x±) equal to the 
classical kink solution, 

h(x±) = m y/6/Ao tanh ^^=M • (2-16) 

As in ( |2.6| ), the kink mass cr( kmk ) can be split into classical, quantum, and counterterm 
contributions: 

a (kink) = a (kink) + a (k ink ) + a (kink) 

By direct evaluation of (|2.7|) , the renormalized part of the classical contribution, which we 
conventionally denote by cr^ mk \ is given by 



a. 



(kink) 
cl 



3 



Ao 

while the explicit form of the counterterm contribution (|2.9| ) is 



4V2^L, (2.18) 



Q./o / .7 \ Jl/2-1 ,,n— 2 
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The quantum piece cr^ 1 ^' can be calculated using (|2.11|) and the spectrum of the operator 
—dj_ + V" . Since we are using dimensional regularization, we can formally rewrite the sum 
in (|2.11|) as an integral, 



a, 



(kink) 
qu 



d" 



(27T) 



dk± 



■(kf + k\ + 2m 2 ) 1 / 2 5'(k ± ) 



[2.20) 



Here d(k±) is the scattering phase shift and the explicit form M of S'(k±) = 



dS(k ± ) 
dk± 



IS 



5'(k±) = -6V2m 



(k]_ + m 2 ) 



(2kl + m 2 )(kl + 2m 2 )' 



(2.21) 



The first two terms in ( [2.20|) come from the two bound states of the small fluctuation 
operator in the kink background, and the last term from the continuous part of the spectrum 
The last term involves a difference between the contributions from the kink and vacuum 
sectors, which do not cancel because of the differing spectral densities in these sectors. 



Substituting ((2.21 ) into ( |2.20|) one arrives at the following expression for cr 



(kink) 



qu 



r (kink) _ ± 



qu 



d n ~ 2 k\\ 
(2tt)"- 2 
dk 



-3\/2m J 
-3V2m 



k 2 



m 



dk ± m 2 (k\ + k\ + 2m 2 fl 2 
2tt (2k ± + m 2 )(k 2 ± + 2m 2 ) 



2tt (k 2 + k\ + 2m 2 y/ 2 (k] 



2m 2 



dk i 



2tt (k 2 + k\ + 2m 2 ) 1 /2 



^2.22) 



Let us analyze Q2.22 ), term by term. The first term that involves (k 2 ) 1 ^ 2 is set to zero 
by dimensional regularization, as it does not contain any external dimensionfull parame- 
ters. The second, third, and fourth terms are finite, and equal 3m/2^/6, — m/v^6, and 
— 3m/v^27r respectively. One way to obtain these results is to formally rewrite each of the 
corresponding (n— 2)-dimensional integrals as a product of an angular integral which is equal 
to 27r n / 2 ~ 1 /T(n/2— 1) and a radial integral which can be shown to contain a (n/2— pole 
that cancels the zero in inverse Gamma function, leading to a finite expression in the limit 
n — ► 2. The last term in ( 2. 22]) can be computed along the same lines, 



3V2 



m 



(2*0 



n-2 



dk | 



3V2- 



m 



2ir {k 2 + k 2 ± + 2m 2 ) 1 / 2 



It obviously produces a pole that cancels the pole in a, 



4tc 



(kink) 



2m 2 

47T/1 2 



n/2-1 



^ n "T(l-n/2) 



^2.23) 



leaving logarithmic dependence 



on the scale \i. This is completely analogous to what happens when one does conventional 
MS regularization in the perturbative sector. 

We are now in position to collect all the terms together, and remove the UV regulator 
by taking the n — >• 2 limit. This leads to the following expression for the kink mass: 



a 



(kink) 



Ay/2 



m(fj,) c 
A 



2y/Q V2tv An 



3^/2 , 2m 2 
+ — — In 



m . 



(2.2A) 
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One can check that this expression is //-independent to 0(A). 

Finally, to make a connection to previous results in the literature, we must establish 
a relation between the MS mass m(/i) and other /i-independent physical quantities. The 
renormalized mass M used in Ref. IJ was defined by the condition that the one-loop effective 
potential has its minima fixed at (0) = ±(6M 2 /A) 1 ^ 2 . This corresponds to the complete 
cancellation of the tadpole diagrams by the counterterm. The relation between M and m(n) 
to one-loop order is given by 

M 2 = m(/j,) 2 + ^- In (2m 2 //j, 2 ) . (2.25) 
Substituting this into Eq. ( |2.20|) gives 



a (kink) = 4^^! + 

A 



M. (2.26) 



which coincides with the well-known result for the kink mass of Ref. • 

The most important lesson that should be learned from this calculation is that it would be 
very difficult to generalize this dimensional continuation approach to the case of an arbitrary 
field configuration. Indeed, the discussion below ( |2.22|) shows that even after a cX . cancels 
the pole in the quantum part of a, taking the n — > 2 limit is still a bit subtle. Namely, it 
involves the cancellation of a zero from the angular integral with a pole from the analytic 
continuation of the radial integral. This procedure would be hard to implement in general, 
as the radial integral for an arbitrary field configuration can only be evaluated numerically. 
Therefore one must consider other ways of regulating a theory, or abandon the zero-point 
energy approach. The former option is chosen in while the latter choice is going to be 
explored below. 



C. Integrating over k± 

One may compute the one-dimensional functional determinant in fl2.10|) directly, without 
converting it to a sum over zero point energies. The benefit, as we have already advertised, 
will be the ability to reduce the computation of the quantum correction to the energy of an 
arbitrary field profile, <r qu + cr c . t . , to a finite integral that can be evaluated numerically. 

It should be clear that the only reason for this procedure not being completely trivial 
is the existence of ultraviolet divergences in the theory. Our aim is the establishment of 
a regulator-independent method, which clearly separates the divergences and allows one to 
use any conventional regularization scheme in dealing with them. 

We will start the analysis of the ultraviolet divergences by considering an expansion of 
the functional determinant in (|2.8| ) in powers of (A — A^)/A^: 



a qu = lim 



1 (-f / A - A"" ' 

2W tr V ~n~ [ A(°) 



^2.27) 



In two dimensions, only the first term in this expansion is divergent. Recall that A^ is just 
an inverse bare propagator, so this divergent term can be written as 



S 



lim [ tr 



{V" - 2m 2 ) (-V 2 + 2m 2 ) 



2\-l 



d 2 k 



2 J (2n) 2 k 2 + 2m 2 ' 
where F is the Fourier transform of (V" — 2m 2 ) taken at zero momentum, 



^2.28) 



dx 



{V" - 2m 2 ) 



— I dx i 
2 



h(x±y 



6m 



[2.29) 



Of course, the integral in ( 2.28|) is logarithmically divergent. However, we are free to 
integrate over the momentum in the direction perpendicular to the worldvolume and then 
impose any convenient cutoff on the remaining integral.^ Subtracting this divergent term 
from Q2.10 ), and adding it back, we obtain the following expression for the finite quantity 
+ o"c.t.: 



cr, 



qu 



(J, 



qu 



+ (T, 



c.t. 



dk\\ 
27 



detj. (k\ 



d 2 + v" 

lim In 

L ^°° det ± (k\ -d\ + 2m 2 
F r^dk\\ 1 



2(kj + 2m 2 ) 1 /2 



^2.30) 



2tt (jfejj + 2m 2 ) 1 /2 



c.t. 



The first integral in this expression is finite, as the divergent part has been subtracted 
from it. The second term in fl2.30| ) contains an integral that appears in the evaluation of the 
tadpole diagram, and needs to be regulated. We indicate this by writing A, which should be 
thought of as some cutoff on worldvolume momentum, at the top of this integral. However we 
can use any regulator we like, as long as it allows unbounded integration over k±. The same 
tadpole diagram also contributes to the mass renormalization in the perturbative sector, and 
consecutively, appears in cr c t . Comparing ( |2.9| ) and ( |2.29|) one can immediately see that 
the UV-divergent terms in (|2.30|) cancel each other, leaving a finite remainder whose value 
depends on the renormalization conditions. 

Thus, the only non-trivial piece of the calculation is the evaluation of a single one- 
dimensional integral. The expression ( |2.30| ) is, therefore, a key result of this section. Note, 
that had the number of dimensions n been greater than two, more subtractions might have 
been needed. The structure of such terms would be quite clear, though, as it is determined 
by the perturbative expansion fl2.27Q . 

In what follows we are going to regulate the UV divergences with dimensional regulariza- 
tion where the number of worldvolume dimensions n—1 is regarded as a variable parameter. 
Since the divergent terms in ( [2.30| ) are isolated, the treatment that is involved is essentially 
the same as appeared when we considered the analytic kink solution. We choose to express 
our result in terms of the physical parameter M defined by ( p. 25 ). One finds simply 



1 r dk\\ 
a ^ + a ^ = 2]^ 



det_ 

lim In 

L ^°° det i 



+ V") 



(k 2 -dl + 2M 2 ) 2(k 2 + 2M 2 y/ 2 



(2.31) 



3 In contrast, the zero-point energy approach does not allow regulator-independent integration over k±, 
an integration over one of the components of kn has already been performed. 



as 
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When the field configuration is a solution of the classical equations of motion, the op- 
erator —dj_ + V" acquires a zero eigenvalue, associated with the spontaneous breaking of 
translational invariance. Consecutively, the 2-dimensional fluctuation operator —d^—dj_+V" 
has continuous spectra extending down to zero. In principle one should go back to the origi- 
nal functional integral and use standard collective coordinate procedures jl^] to separate the 
translational mode. However, the zero mode of —&]_ + V" merely causes the two-dimensional 
functional determinant appearing in fl2.31p to have an integrable logarithmic singularity at 
fcn = 0. Therefore, we do not need to modify ( [2.31| ) for the computation of the one-loop 
kink mass. 

To proceed further, one needs to evaluate the ratio of one-dimensional functional deter- 
minants appearing in ( p.31|) . Fortunately, we can use the result of [|13|] where this ratio is 
directly related to solutions of associated differential equations. (We review this construction 
in the appendix.) Specifically, this ratio of one-dimensional determinants can be written as 



detg (fcjj -dj + V") ip(L) 
det ± (kj-d\ + 2M 2 ) ~ V (0 W 

where the t/>'s are the solutions to the associated differential equations, 



(2.32) 



*f " dl + V"] 4(x ± ) = 0, 



(2.33) 



and 



k\ -dl + 2M 2 ] ^ {0) (x ± ) 



(2.34) 



with initial conditions 



<K-L)=0; ip(—L)' 



(2.35) 



and 



^ (0) (-L)=0; ^°>(-L)' = l. 
The free solution ip( \x±) is trivial to find. It is given by 



(2.36) 



^°\x ± ) 



sinh 



{kf l +2M 2 ) 1 / 2 {L + x ± ) 



(fcjf + 2M 2 )V2 



(2.37) 



Neglecting terms which are exponentially small as L 
to (|3D is 



oo, the contribution of this solution 



-(k 2 +2M 2 ) 1 / 2 2L 

e II 

2(A; 2 + 2M 2 )V2- 



(2.38) 



It is desirable to extract explicitly the large L asymptotic behavior from ip{x±) as well. 
This prompts the natural substitution 



ip(x±) = f(x ± )e~^ 



(fe2+2JW 2 ) 1 /2( a; _ L+ I ; ) 



(2.39) 
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and converts the ordinary differential equation (|2.33|) to the corresponding equation for 



~-d\ + 2{k\ + 2M 2 ) l ' 2 d L + (V" - 2M 2 )] f(x ± ) = 0, 
with f(x±) being subject to the same initial conditions, 

f(-L) = 0; f'(-L) = 1. 



(2.40) 



(2.41) 



We can now present the result for the quantum correction to the mass of an arbitrary 
kink-like field configuration in two-dimensional 4 theory in its final form, 



1 r dk 



Vqu + Cc.t. 



(27T) 



lim In 

L^oo 



2{k\ + 2M 2 ) 1/2 f{L) 



2{k 2 + 2M 2 )V2 



(2.42) 



As a check of our method, we have used ( 2.40|) -( p.42|) to recalculate the value of cx qu + cr, 
for the kink configuration, where the exact result is given by (see Eq. fl2.26|) ) 



at. 



(kink) , (kink) 
°qu ' a c.t. 



2VQ V2 



71 



M. 



(2.43) 



Using the canned routines of Mathematica to integrate Eq's. (2.4U)-( 2.42p , one may easily 
reproduce this result with an accuracy better than one part in 10 5 . 



III. THREE-DIMENSIONAL SCALAR ELECTRODYNAMICS 

We are now ready to tackle the gauge field theory case. For the rest of the paper we 
will be developing a method for computing one-loop energy densities (which will sometimes 
be called the surface or interface tension) of domain wall-type Higgs field configurations in 
three-dimensional scalar electrodynamics. The parameters of this theory will be tuned to 
sit at the first order phase transition where, at one-loop level, there is a degeneracy of the 
ground state, and consequently the existence of stable domain wall-type configurations. 

As we will see, the method that we applied in the previous section can be successfully 
generalized to scalar QED. However, there are two new problems which we did not en- 
counter in the scalar theory. One is the existence of quadratic couplings between different 
fluctuating fields in the Lagrangian. This requires generalization of the result for a one- 
dimensional functional determinant. Another issue is the absence of tree- level degeneracy 
of the effective potential which will cause the quantum and classical parts of the surface 
tension to individually diverge with increasing box size. We will show how these problems 
can be resolved. 



A. Formulation of the problem 

We will start by formulating the theory and discussing the relations between the pa- 
rameters that are required to tune the theory to its phase transition. We will see that the 
parameter that determines the validity of perturbation theory is the ratio of scalar and gauge 
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couplings. In the regime where this parameter is small, it is in fact possible to analytically 
compute the leading behavior of the interface tension. We will later compare the results of 
this analysis with the results of our full one-loop computations. 

We consider the following Euclidean action for three-dimensional scalar electrodynamics, 



S = d 6 x 



^FijFij + ^(A$) f (A$) + V(Q) + £ g . f . + £ gh + £ c .t. 



(3.1) 



where A = d\ + igA{ denotes the usual covariant derivative with respect to the abelian 
gauge field A iy and £ g .f. and £ g h denote gauge fixing and ghost contributions to the La- 
grangian, respectively, which will be specified later. We will eventually be using dimensional 
regularization, so the counterterm Lagrangian £ ct would in general contain (n— 3) _1 poles 
where n is the analytically continued number of dimensions. In three-dimensional scalar 
QED, however, no poles appear at one-loop order. Hence, we may drop explicit mention of 
counterterms in all subsequent calculations. 

We choose the simple quartic potential for the complex Higgs field $, 

vm = ^-m 2 + ±m\ (3.2) 

with a positive coefficient m 2 in front of the quadratic term. 

Let us briefly sketch the relations between the parameters of the theory that follow from 
the condition of a doubly degenerate ground state. Recall that because of the U(l) gauge 
symmetry, the effective potential V e s may be chosen to be a function of a real expectation 
value of the Higgs field. The crucial observation is that for the theory at hand, the effective 
potential acquires a (— g 3 $ 3 ) term which comes from the summation of all one-loop gauge 
field diagrams with an arbitrary number of insertions of g 2 $ 2 [14]. The existence of two de- 
generate minima of the effective potential hence implies that the effective potential acquires 
the form of a double well, 

2 

777 

V e ^— v) 2 , (3.3) 

where v is the expectation value of $ in the Higgs phase. This immediately implies the 
following characteristic relations between the parameters: 

v ~ g 3 /X, (3.4) 
m 2 ~ g e /X. 

Given ( |3.4|) , one can verify that the parameter controlling the reliability of perturbation 
theory in both the symmetric and Higgs phases is X/g 2 |[L4| . For example, adding an extra 
loop in a diagram in the Higgs phase due to the addition of a gauge field propagator leads to 
an additional factor of order g 2 J d 3 k(k 2 + Mjy)~ 2 ~ g 2 /M w ~ X/g 2 , where M w = gv is the 
mass of a gauge boson. Thus, perturbation theory is reliable for small X/g 2 , even though 
the existence of a Higgs phase is not apparent at tree level, and is caused only by quantum 
corrections. 

As in the previous section, we put the system into a three-dimensional box of half-length 
L and transverse area V». The computation of the one-loop surface tension amounts to 
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the evaluation of the one- loop effective action for a given field profile h = h(x±), which 
we may choose to be real and positive. We will restrict our attention to domain wall-type 
configurations that obey boundary conditions similar to ( |2.2| ): 



h{-L) = 0, 
h(L) = v. 



(3.5) 
(3.6) 



We decompose the complex Higgs field $ into the background field h plus two real fluctuating 
fields, 



$ = h + x + i<f>, 

and attribute gauge variations to the fluctuating fields x an d 4>- 
The effective action is a sum of classical and one-loop terms: 



(3.7) 



T[h] = T d [h] + T qu [hl 



where 



r c i[/i] = V\\ J dx< 



I 777 ^ \ 

-(d ± h( x± )) 2 + — h(x ± y + ^h( x± y 



(3.8) 



(3.9) 



To get the quantum part, we need to expand the action (|3.1|) around the background field 
configuration, keeping terms that are quadratic in fluctuating fields. In doing so, one notices 
that a term of the form Aidi<p appears. In the presence of a spatially varying background 
field, it is not possible to eliminate entirely this cross-coupling term. But one may partly 
eliminate this term by choosing a modification of gauge with the following gauge-fixing 
term, 



«9^ + £f Im($ 2 ) 



The part of the action quadratic in fluctuating fields then reads 

5( 2 ) = J d 3 x (C x + C A + + C A(t> + £ gh ) , 



where 



£ = ^0) 2 4 



X 



-Ai {-5 l3 d 2 - (1 - C X Wj + g 2 h(x ± ) 2 5 tl ) A, 



m 



+ 



A £ 2 

h -g 

12 T 



h(x ± ) 2 cf> 2 



A4> 



-2g(d±h(x ± ))A. 



C^ = f 1 {-d 2 + ig 2 h{x i _) 2 )r 1 . 



(3.10) 

(3.11) 

(3.12) 
(3.13) 

(3.14) 

(3.15) 
(3.16) 
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It is particularly convenient to choose £ = 1, because with this gauge choice the terms 
containing (VA) 2 disappear. Notice that only the A± component of the gauge field is coupled 
to the field, thanks to the gauge choice ( |3.10| ). 

Performing the functional integral gives a product of functional determinants, as always. 
One can verify that with the gauge choice that we have made, the contribution of the ghosts 
completely cancels the contribution of the two A\\ components of the gauge field. 

The expression for r qu is thus composed out of contributions from the Higgs boson \ 
and the coupled fields A± and </>: 

r qu [h] = r x [h] + r M [h], (3.17) 

F x [h] = In det ( ^) , (3.18) 




T M [h] = In det , (3.19) 



where the fluctuation operators can be read off directly from (|3.12j )- (|3.16| ). Before we write 
them explicitly, let us introduce notation that will simplify the subsequent formulae, by 
defining 

M H {hf = m 2 + h 2 , (3.20) 
M w {h) 2 = g 2 h 2 , (3.21) 

M A ^h) 2 = m 2 + (^ + g 2 )h 2 . (3.22) 
o 

With this notation, the operator that determines quadratic fluctuation of x is 

A x = -V 2 + M H (h) 2 , (3.23) 

where according to ( |3.20| ), Mh{K) 2 = m 2 + ^h(x±) 2 . The fluctuation operator that couples 
A± and is denoted by Aa</>- It is a linear operator on TCl®'Hl 1 and can be represented 
by the 2x2 matrix: 



-V 2 + M w (h) 2 -2g(d ± h(x ± )) \ 

-2g(d ± h(x ± )) -V 2 + M^(^) 2 I" {6 ^> 



In Eq. ( |3.19| ), A^i denotes the corresponding operator in the symmetric vacuum sector .0 

We can rewrite the expressions (|3.18|) and fl3.19| ) using the fact that the dependence 
on worldvolume coordinates of these operators is trivial, and thus they can be partially 
diagonalized by a Fourier transform.^ That is, 



4 We could equally well have chosen to be the fluctuation operator in the Higgs phase, since the two 
ground states are degenerate 

5 This is what we also did in section 2, by rewriting er qu as an integral over the worldvolume momentum. 
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2 
2 



dkf 1 
(2vr) n - 1 



detj_ A v 
In x 



det ± 



(27T) 



n-1 



In 



det i A 



A(j> 



det i A 



(0) ' 



where 



and 



A 



A x = k \\ - d l + M H{hf 



-dl + kl + M w {h) 2 -2g(d±h(x±)) 
-2g(d ± h(x ± )) -d\ + k 2 + M M (h¥ 



(3.25) 
(3.26) 

(3.27) 
(3.28) 



As we saw in the previous section, the divergent integrals over the worldvolume momen- 
tum in (|3.25|) and ( |3.26|) can be regularized using any convenient regulator. We just need to 
subtract the divergent pieces and add them back in a regularized form. We find it convenient 
to use dimensional regularization. 

One can compute the effective potential by taking the field configuration to be coordinate- 
independent and going to the limit L,V\\ — > oo. Then the operators A x and Aac/> can be 
further diagonalized, giving the following expression for the effective potential: 



d n k 



In 



'k 2 +M H ($y 
k 2 +m 2 



In 



'k 2 +M w ($y 
k 2 



In 



'k 2 +M A ^f 
k 2 +m 2 



. (3.29) 



The integrals can be easily evaluated using dimensional regularization. The resulting 
effective potential (which has also been computed in many references, see e.g. Refs. [fT , 10 , lT| ) 
reads 



(3.30) 



As we asserted earlier, a cubic g 3 $ 3 = Mw(^) 3 term has appeared, leading to the 
relations (|3.4|) for the parameters of the theory. 



B. Computation of the surface tension 

The energy density of a domain wall-type configuration is a sum of the classical and 
quantum parts, as defined in (|3.8| ). Both classical and quantum terms immediately follow 
from ( |3.9| ) and ( |3.17 ), by dividing the corresponding expressions by the transverse area 
V\\. Note, that a c \ is now linearly divergent with increasing L, as there is no ground state 
degeneracy in the tree-level effective potential. This means that the quantum part of a must 
diverge with increasing L as well, so that the total a is finite as L — > oo. Hence, we must 
take the L —>■ oo limit at the end of the calculation, after all terms are collected together. 
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The quantum part, <r qu , can be further decomposed analogously to (|3.17|) with corre- 
sponding terms called a x and cx^, respectively. 

Let us concentrate on the calculation of a x first. We would like to apply the same 
technique for computing the ratio of determinants as used in the previous section. However 
the situation is now different, since v = lim h(x±) is no longer the minimum of the tree- 

level effective potential, and thus the ratio of determinants diverges with L instead of being 
asymptotically constant. In fact, it must diverge exponentially, so that after taking the 
logarithm, the leading divergence is linear in L. 

Naturally, we would like to extract this asymptotic behavior, leaving to the computer 
only a calculation of a term that is finite in the limit L —>■ oo. For that purpose, we impose 
a restriction on our field configuration. We require it to be constant outside the region 
—a < x± < a; so that h(x±) = for x± < —a and h(x±) = v for x± > a. In other words, 
we require the domain wall field profile to have a "thickness" no larger than 2a. Eventually, 
a will be taken sufficiently large so that this condition has negligible effect. 

Thus there now are two length parameters that we have introduced: a, being the half- 
length of the window where the background field is allowed to vary, and L, the half-length 
of the box in the direction perpendicular to the worldvolume. The quantum part of the 
surface tension, as well as a x and oa<$> individually, are going to be asymptotically constant 
in a, but linearly divergent in L as we send L to infinity. 

The ratio of determinants in a x can be written in terms of the solution of associated 
ordinary differential equations, just as in Eq. (|2.32|) . This leads to the following expression 
for o x : 

1 f dn " lfe J In ( \ 



a * = 2j (^p ln (w' (3 - 31) 

where ip(x±) is the solution of the differential equation 

'-d{ + kl + M H (h) 2 ] rf;(x ± ) = 0, (3.32) 
with the boundary conditions 

?/>(-£)= 0; ij){-L)' = l. (3.33) 
The corresponding solution in the vacuum sector, ip^(x±), can be taken from ( |2.38| ): 

{k^+m 2 ) 1 ' 2 2L 

Note that Eq. ( |3.31| ) must be rendered UV-finite by the separation and subsequent di- 
mensional regularization of the UV-divergent piece. At this point we only note that the 
regularization proceeds according to the procedure established in section 2. We will return 
to the actual treatment of ultraviolet divergences in a x later. 



To solve the differential equation (|3.32|) we divide the whole interval into three parts. 



Namely, region I: x± < —a, region II: —a < x± < a, and region III: x± > a. We can 



immediately write the general solution of (|3.32| ) in regions I and III, since h(x±) restricted 
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to these intervals is constant. Thus, we only need to integrate Eq. ( |3.32| ) in region II, and 
then match the solution at the intersections x± = —a and x± = a. 

Omitting terms that fall off exponentially with increasing L, the solution in region I is 
given by 



1 



(fc 2 +m 2 ) 1/2 (L+x±) 



2(fcg + m 2 ) 1 / 2 

The solution in region II can be written in the following form: 

(fc 2 +m 2 ) 1/2 (L+x±) 



x± G [-L, a] 



(3.35) 



V>(z_l) 



2(fc 2 + m 2 )V2 



x± G [—a, a] 



(3.36) 



where f^ n \x±) satisfies the equation 



-di-2(ki + m 2 yi 2 d± + h(x 1 _y 



f {H) (x ± ) = 0, x ± E[-a,a] 



(3.37) 



with the boundary conditions 

/(")(-«) =1; /W(-o)' = 0. (3.38) 
The general solution in region III is 

^( x± ) = j(m) e (fc 2 W)V 2 L+(fc 2 +Af H H 2 )V2 a . ± + ^ (m)e ( fc 2 +m 2 ) l/ 2L _ (fc 2 +MHM 2 ) l/2 x± ^ ^ ^ 

(3.39) 

Matching this solution with ( |3.36| ) by requiring the continuity of ip(x±) and its first deriva- 
tive, we obtain the following expression for 



f 



(in) 



a[(fc 2 +m 2 ) 1 / 2 -(fc 2 +M w (i)) 2 ) 1 /2] 



f {II) {a)' 



: I l "ill + f( H )( a ) 

2{k 2 + m 2 ) 1 / 2 \ {k 2 +M H {v) 2 y/ 2 J y) 



(k 2 + m 2 ) 1 / 2 

1 + (k 2 +M H ( v yy/ 2 



(3.40) 



It is f( IH ^ that determines the value of a x , as the second term in (|3.39|) is exponentially 
suppressed in L. One can therefore take the first term in ( |3.39| ) with f^ IH ^ given by ( |3.40| ) 
and substitute this, together with ip^°\L) from ( |3.34j) into the expression ( |3.31| ) for a x . 
Before we do that, however, note that for the constant background field h(x±) = v, ip{L) 
can be easily found analytically. This allows one to compute the part of the effective action 
T x [h = v] at finite L: 



T x [v] _ 1 f d n ~ x k 



2 J (2n) n - 1 



2L 



(k 2 + M H (vyyi 2 -(k 2 + m 2 y/ 2 



1 fk 2 + M H (vf 



+ - In 

2 



kn + m 2 



(3.41) 



The term that multiplies L in this expression comes from the exponent in ( |3.34j ) and the 
corresponding exponent in ip(L), while the ratio inside the log in the last term involves the 
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coefficients that multiply the respective exponents. One can obtain part of the effective 
potential [the M#($) 3 term in Eq. (|3.30|) 1 by dividing the expression above by the volume of 
the spacetime box 2V\\L and taking the limit L — > oo. The last term in ( |3.41| ) is a boundary 
effect reflecting our choice of Dirichlet boundary conditions. It makes no contribution to the 
effective potential in the L — ► oo limit. 

We may require the parameters of the theory formulated in a box with finite length L 
to be adjusted so that the effective action evaluated at finite L vanishes at v. Of course, as 
we send L to infinity, these parameters will flow to the corresponding values for the theory 
in infinite spacetime. This condition will prove useful momentarily. For one can notice that 



the part of cr x Vi\ that is linear in L is simply half of the first term in ( 3.41 ). Thus, one can 
separate a x as 



0\ 



2VFi 



(3.42) 



where 



d n - x k 



(2n] 



n-l 



In 



2 j(m) 



(H + m 2 )V2 



-iln 

4 



k\ + M H (v) 



hi 



(3.43) 



is finite as L — > oo and T x [v] is given by Eq. ( |3.41| ). 

Recall that T x [v] is a part of the effective action evaluated at the constant field configu- 
ration h(x±) = v. One can anticipate that the same story will repeat for the calculation of 
o"A0, and r^fi;] will be produced in a similar fashion. Thus, after we collect all parts of a 
together, the coefficient in front of L in the resulting expression will be a multiple of T[v], 
and therefore the divergence in L will vanish, leaving an expression for a that is finite in the 
limit L — > oo. 

Having anticipated the form of the outcome, let us now compute cr^. The relevant 
differential operators in (|3.19|) are now given by 2 x 2 matrices, so we must find a way to 
generalize the correspondence between determinants and solutions of associated differential 
equations. This is done in the appendix, where we show that 



det A M _ ^\L)^\L)-^>{L)^\L) 



',(2)/ 



det A?) 



<pP(L)<pP(L) 



(3.44) 



where ift^(x±) and ip^(x±) are solutions of the same system of differential equations, 



A 



AS 



4\ XA _) 



0. 



1, 2, 



(3.45) 



with vanishing initial conditions, 




1, 2, 



(3.46) 



but with differing initial derivatives, 
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1 





4 2), ( 





1 



(3.47) 



In other words, tp^ and ip^ are the two linearly independent solutions of A^^W 
which vanish at x± = —L but which have orthogonal initial "velocities" . 

The corresponding vacuum sector solutions are denoted by cp®. 

Hence, <JA<t> can be written analogously to ( |3.31| ) as 



dk 



n— 1 



$\L)W\L)-W\L)W>{L) 



,(2) 



0-h 



',(2), 



111 



(3.48) 



The further computation of proceeds similarly to that of a x in almost all respects. We 
again divide the interval [-L, L] into three regions, and use appropriate changes of variables 
to extract the asymptotic behavior of the solutions. The main results are presented below. 

As in the case of a x , ga^ can be decomposed analogously to ( |3.42|) : 



PA<t> — ® a<$> + 



A<f> 



2V\\ 



(3.49) 



where 



&A4> 



(2k) 



n-1 



In 4 



Al,III) f {2,111) f {2,111) f {1,111) ' 
Jl J 2 ~J1 J 2 

(kj + m 2 ) 1 / 2 k\\ 



-iln 

4 



'(k 2 +M w (v) 2 ) (k 2 +M M (v) 
k 2 (k\ + m 2 ) 



(3.50) 



and 



(2k) 71 - 1 



2L{(k 2 + M w (vf) 1/2 -h + ( k l 



(k 2 + M w (vf ) (k 2 + M M (v) 



1 . 
+-ln 
2 



k 2 (k 2 + m 2 




v) 2 ) 1 ' 2 -(k 2 + m 2 fl 2 } 
(3.51) 



where /j is determined analogously to f^ 111 ^ in ( |3.40| ). The subscript i in the former 
labels the component of a vector, while the superscript j distinguishes the two sets of 
boundary conditions. In other words, fi' 11 comes from the if)!*? (x±). 
Explicit expressions for fj; 1 ' 11 are: 



fi 



{i,ni) 



{i,ni) 



o a[k\ r {kl+M w {v¥Y'*} 

2k H 



3 o[fc||-(fcJ+M^(«)»)Va] 

2k n 



fi^'ia) 



(k 2 + M w ( v ) 2 y/< 



+fi 1 ' I %) 



f, 



{hay 



(k 2 + M A ^(v) 2 yr< 



+f: 



{hii) 



(k 2 + M w (v) 2 y/ 2 



(k 2 + M M (v) 2 y/ 2 



(3.52) 



(3.53) 
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where /^'^(xj.) and f2 1 ' II \x±) are the solutions of the following system of ODE's, 

A (M/) \ 

subject to the boundary conditions 



-d\ - 2A;j,9 ± + M w {h) 2 -2g(d±h) 

-2g(d ± h) -d 2 ± - 2k\\d ± + M A<p (h) 2 



0, 



(3.54) 



f?' H \-a) 
ff JI \-a) 







f, 



-a) 








(3.55) 



Corresponding expressions for /, 



(2,111) 



are very similar, as one might expect. The only 



difference is caused by the fact that tp^(x±) and il)^ 2 \x^_) have different exponents in the 
region I. One finds, 



ft 



(2,111) 



f: 



(2,m) e 



o a[(kl+ m ^-(kl+M w (v)^] 

2(k 2 + m 2 ) 1 / 2 

«[(fc5+m 2 )'/2-(fcJ + M^(„)2)V2] 



A (2 ' 7/)/ (a) 



(H+M w (t;)2)i/2 



f/f'^V) I 1- 



(ifej + m 2 ) 1 / 2 

\k 2 +M w (vyy/ 2 



(3.56) 



2fcn 



/ 2 (V/) '(a) 



(A; 2 +M A0 (t;) 2 )V2 



(k 2 + m 2 ) 



2U/2 



(A;g+M^(t;) 2 )^ 

(3.57) 



where /f' n '(xi) and /^''"''Oex) are the solutions of the following system of ODE's, 



e{2,II) 



-d 2 ± -2(kj+m 2 ) 1 / 2 d ± -m 2 +M H (h) 2 -2g(d±h) 

-2g(d ± h) -dl-2(k 2 +m 2 ) l l 2 d ± -m 2 +M A4> (h) 2 



ft 
f: 



(2JI) 



(2,H) 



o, 



subject to the boundary conditions 

( fi 2 ' H \-a)\ . ' 
\tf n \-a)) " 







fi 2JI) \-a) 







(3.58) 



(3.59) 



Now one can collect together a c \, a x and cr^, and in particular combine the terms 
which add up to make r[t;]/2Vj|, and set it to zero. At this point we need to recall our 
prescription for regularization, which was introduced in section 2. Namely, we need to 
subtract the divergent pieces from a x and a A ^>, and then add them back, with the divergent 
loop diagrams replaced by their regularized values. 

Naively, the divergent term in a x , for example, is given by 



4««v> = ± G tr 
x 4 



(V 2 + m 



2\-l 



where 



G = J dx±h 2 (x±) 



(3.60) 



(3.61) 
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is defined analogously to F in Eq. ( |2.29| ), and the trace involves integration over k» and 
summing over discrete x± modes that vanish at ±L. However, when extracting r[v]/2V|i, 
we should also add and subtract the corresponding UV-divergent term which is just 



A r 2 

-Lv 2 tr 
4 



(V 2 + m 



2\-l 



A 
4 



dx±v tr 



(V 2 + m 



2\-l 



(3.62) 



Therefore, the UV-subtraction in a x should actually involve 4 div ) as in Eq. flOOD but with 
G replaced by 



G 



dx 



h 2 (x ± ) -9(x ± )v' 



(3.63) 



Here 6(x±) is the step function which is zero for x± < and one for x± > 0. Note that G 
is finite in the limit L — > oo. 

The treatment of ultraviolet divergences in proceeds in a similar manner. 

Taking the L —>■ oo limit we arrive at an expression for a that only involves finite integrals 
that can be computed numerically: 



a = lim [c^i + cr x + a A( p] , 

L^oo 



(3.64) 



where 



0". 



if '[ 771 r 

<fe_L < -(d±h(x±)) 2 + — [/i(x ± ) 2 



+ 4! 



rf 2 fc|l 



In 



2 j(///) 



(k 2 +m 2 y/ 2 



- - In 

4 



fc 2 + m 2 



/i(x ± ) 4 - #(x ± )t/ 



G 



4 (fcjj+m 2 ) 1 / 2 



A G 
+-— m, 

4 47T 



/■ d 2 k\\ 




In ^ 


1 (2tt) 2 















,(l,m) , (2,7/7) , (2,///) , (1,7//) 1 
I J 1 J2 Jl J2 

{k 2 + m 2 ) l l 2 k\\ 
' (k 2 +M w (v) 2 )(k 2 +M w (v 2 ))\ g 2 G ,A g 2 



(3.65) 



(3.66) 



G 



A; 2 (A; 2 + m 2 



.A g 2 .G 

+ ( h — ) — m. 

42 2 ; 4vr 



2 fey v 12 2 7 (A; 2 + m 2 )V2 

(3.67) 



Note that the "classical" term ( |3.65| ) is also finite as L — > oo, as it is defined by the same 
subtraction that was needed to isolate r cl [-u], 



0"cl = Ccl 



2Vn 



(3.68) 



The terms outside the integrals in a x and are the values of UV-divergent terms of the 
type ( p. 60 ), computed in the limit of infinite L by dimensional continuation. 
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C. Numerical evaluation 



It is now possible, by using Eq's. ( |3.64] )-( |3.67D , to compute the full one-loop surface 
tension for an arbitrary domain wall-type field configuration. One just needs to solve the 
relevant differential equations ( |3.37|) , (|3.54|) , and ( |3.58|) , and then perform a final integration 



over the worldvolume momentum k\\. The latter step merely involves a one-dimensional 
integration over the radial component of kn, as the integrands in a x and depend only 
on 

We have written a simple Mathematica program to implement these steps. It allowed 
computation of a with a relative accuracy of 0.2% or better .Q As described below, we tested 
the numerical results against perturbative calculations when X/g 2 << 1, and obtained good 
agreement. 

In principle, one may calculate the one-loop energy of the true domain wall configuration 
by minimizing ( J3.K4 ) with respect to the domain wall profile. This we have not bothered to 
do. 



D. Comparison with perturbative X/g 2 expansion 

As noted earlier, for our specific example of the Abelian Higgs model, the ratio of scalar to 
gauge couplings, X/g 2 , functions as the loop expansion parameter controlling the reliability 
of perturbation theory at the phase transition. For asymptotically small values of X/g 2 one 
may analytically compute the leading and first subleading correction to the surface tension 
a. One way to do this is to start with the derivative expansion of the effective action: 



d x 



V cS (h) + ±Z(h)(d,h) 2 + - 



(3.69) 



where • ■ • stands for the terms with four or more derivatives. We denote the first two terms 
in this expansion by r< 2 >. 

To leading order in rj = X/g 2 one may take the wavefunction renormalization factor Z{h) 
in the effective action (|3.69|) to equal one. The reason for this is that the one- loop corrections 



to the classical kinetic term in ( gjfjg ) behave like g 4 h 2 M w (h)~ 3 (dh) 2 and X 2 h 2 M H (h)- 3 (dh) 2 , 
and hence their contribution to a is suppressed in comparison with the leading term by 
a factor of g 2 /Mw{v) ~ r\ and X/Mh(v) ~ rfl 2 , respectively. (Recall, from (|3.4j ), that 
Mw(v) — gv ~ g 4 / X while Mh(v) ~ m ~ g 3 / VX.) Hence, the leading term in the expansion 
of a is completely determined by the leading term in the expansion of the one-loop effective 
potential in powers of rj. Recall that for a scalar theory with a canonical kinetic term and 
any effective potential with two degenerate minima at $ = and $ = v, normalized so 
that V e s(0) = V e g(v) = 0, the domain wall field configuration is a solution of the classical 
equations of motion whose conserved "energy" E = h(d±h) 2 —V e s(h) vanishes. Consequently, 



6 Higher accuracy is undoubtedly feasible given more care t han we used in the numerical implementation. 
This is not completely trivial however, as Eq's. ( |3.54 ) and ( 3.5§| ) can be quite stiff due to the presence of 



two different length scales set by the mass parameters Mw(h) and M^^/i). 
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the surface tension, neglecting the wavefunction renormalization and all higher derivative 
terms, is given by 



a 



dx ± (d ± h) 2 = ( V dh j2V eS (h). (3.70) 
Jo v 



To obtain the leading behavior of the integral ( |3.70|) , one may expand the one- loop 



potential ( |3.30| ) for small 77, obtaining just the classic double well form ( |3.3| ). Inserting this 
into Eq. (3.70) gives 

a = *® (1 + 0(77)) , (3.71) 



a (o) = ™ T dh h (v _ h ) = \ mv > = _A. ^-3/2. (3.72) 



with 



v Jo 6 v^tt 

[The final form uses the leading behavior of v 2 = ^~ 3 / 2 777, (1 + 0(rj)).\ To get the com- 
plete next-to- leading term, it is necessary to calculate both the two-loop effective potential 
and the one-loop correction to the wavefunction renormalization, i.e., Z(h) — 1. One can not 
go beyond the next-to-leading order term using the derivative expansion, since the Higgs 
field contributions to the higher order terms in the expansion ( |3.69| ) give corrections to the 
surface tension that all scale like a^rf/ 2 . To see this, one may look, for example, at the 
contributions to the two- and four-derivative terms that behave like X 2 h 2 Mi{(h)~ 3 (dh) 2 and 
A 2 M#(/i)~ 5 (<9/i) 4 , respectively. The first term represents a Z-factor correction and scales 
like a^r/ 3 / 2 . The second term is further suppressed by a factor of [IMh(v)]~ 2 where / is the 
characteristic thickness of domain wall. But this is 0(Mh(v)~ 1 ) and so this four derivative 
contribution also scales like a^rf' 2 . It is easy to see that further higher order derivative 
terms also scale like a^rfl 2 . 

More generally, the expansion of a in powers of 77 breaks down because of infrared 
problems M. This can be seen by looking at the contribution of the gauge field to the four 
derivative term computed to one-loop order. This term behaves like g 4 M\y(h)~ 5 (dh) 4: and 
by naive power counting gives a contribution to a that scales as a^'rf. However, since h(x±) 
vanishes exponentially in the symmetric phase (i.e., as x± — > —00), and Mw(h) = gh, the 
integral over x± in the effective action diverges. 

In fact, there is another, more significant source of infrared divergences. One can see 
this by expanding the effective potential to higher orders in 77. One finds that the integral 
which enters the 0(rf/ 2 ) term is logarithmically divergent in the infrared. In particular, this 
implies that the next-to-next-to-leading order correction to a is 0(r] 3 ^ 2 log 7]). 

Nevertheless for small X/g 2 , one might hope that the one-loop value of a could be well 
approximated by the contributions from the effective potential plus the first wavefunction 
renormalization correction, while neglecting all higher order terms in the derivative expan- 
sion. It was argued in Ref. that this is indeed the case in electroweak theory, at least 
for intermediate values of the Higgs masses. Unfortunately the infrared divergences of the 
individual terms in the various expansions used in [[J complicate the accurate computation 
of a. In contrast, the technique that we use in this paper does not have infrared problems, 
as it effectively resums all one-loop corrections to the surface tension, allowing its accurate 
computation for any value of X/g 2 . 
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In what follows we show that the inclusion of the wavefunction renormalization gives 
results that are very close to the full one-loop value of a given by ( 3.64|) , as long as \/g 2 is 
less than 1. The comparison with analytic results also provides a nice consistency check of 
Eq. flpl). 

To make this comparison, it is necessary to calculate the wavefunction renormalization 
in the particular gauge ( [3.101) . We use the convenient method of Ref. [15| to determine the 
first two terms in the derivative expansion of the effective action ( |3.69j ), which we denote by 
I^ 2 ). According to Ref. ||15||, the one-loop contribution to can be written as 



5T 



(2) 

one— loop 



ST^ + 5T 



(2) 
A<j>i 



where 



d X 



d 3 k 



n — — A 

Ai 0) + 3 



(3.73) 



(3.74) 



and 



5T 



(2) 



d x 



d 3 k 



tr ln xw + y tr (^ A i 



(3.75) 



In the equations above, the trace runs over matrix indices and the hat on top of the operators 
A x and [defined in Eq's. ( |3.23| ) and ( f3.24| )] denotes a Fourier transform with respect 
to k. For example, 



A x = k 2 + M H {h) 2 . 
The wavefunction renormalization 

Z{h) = 1 + 5Z X + bZ A4> 

receives contributions both from 5T^ and ST^l, denoted by 5Z X and 5Z/ 
The scalar contribution can be easily evaluated: 



\ 2 h 2 



5Z^ 



d 3 k 



k 2 



X 2 h 2 



(2vr) 3 (k 2 + M H (h) 2 Y M H (h) 3 192vr 



(3.76) 

(3.77) 
respectively. 

(3.78) 



To calculate ^r^j, one should notice that contributions with two derivatives come from both 
terms in Eq. ( |3.75| ), as the fluctuation operator ( |3.24|) already contains the derivative term 
g(dh). A computation which is analogous to the one done in f9j leads to the result, 



9 A h 2 



5Z 



A<t> 



+ 



M w (h) + M A ^h)n \M w (hf M A<p (h) 



48tt 



(3.79) 



One can easily check that SZa^ ~ rj and 5Z X ~ rj ' , in accord with our earlier discussion. 

The expression for the surface tension which takes into account the Z-factor is a simple 
generalization of ( |3.70| ). The vanishing "energy" now implies that j)Z(h)(d ± h) 2 = V eS {h), 
and expression ( |3.70|) becomes 
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FIG. 1. Co mpari son of the full one-loop surface tension (3.64) [solid line], the second order derivative 
expansion result ( [3.8C ) computed wi th th e full one- loop effective potential and wavefunction renormalization 
[long dashes], the "classical" result (3.7C) computed with the full o ne-loop effective potential [short dashes], 
and the first two terms of the small r\ asymptotic expansion (3.84) [dotted line]. 



a 



dx ± Z(h)(d ± h) 2 = [ V dh ^2Z(h)V eS (h). 



(3.80) 



We can compute the first two terms in the expansion of this zT-factor corrected a by plugging 
in expressions ( |3.77| )-( p.79| ) for Z(h), and ( |3.30| ) for V c s(h), rescaling the field as h = vh to 
extract a factor of v, and expanding the integrand in powers of rj. The leading term in the 
expansion is just a^°\ To get the subleading term it is necessary to use values of g and v 
accurate to next-to-leading order: 



g = S^vW 74 ™ 172 



32 



v 



48 n i/4 



1 



" 3/4 m 1/2 (1-^1, 



0(rf' 2 ) 



One finds that the first subleading term is given by 



r(l) 



m 



8^' 



-1/2 



dh (h - 1) (17+ 18/i) 



23 m 2 



-if 



-1/2 



71 JO 16^^ 

Hence, the zT-factor corrected surface tension ( p.80j ) has the small rj expansion 

23 



a 



a 



(0) 



1 



32 



7] + 0{rf 12 log 77) 



(3.81) 
(3.82) 

(3.83) 

(3.84) 



with given in Eq. (gTg) . 

The result ( |3.84j ) is not the complete next-to-leading order surface tension, as two-loop 
corrections to the effective potential also contribute to the subleading 0(rj) term. Including 
these two-loop contributions is straightforward, but will not be examined here. 

It is interesting to compare the different approximations to the surface tension given by 
(i) the full one-loop result (|3.64j ), (it) the expression (|3.80|) computed with the complete 
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one-loop effective potential ( |3.30|) and wavefunction renormalization factor ( |3.77|) , (Hi) the 



expression (|3.70| ) computed with the complete one- loop effective potential ( |3.30|) , and (iv) 



the first two terms of the small 77 asymptotic expansion ( |3.84| ). All these expressions omit 
two (and higher) loop contributions to the effective potential, while (ii) also omits four 
derivative (and higher) contributions to the one-loop effective action, and (Hi) omits all 
one-loop derivative corrections to the effective action. This comparison is shown in Fig. |l|. 
The full one-loop result ( |3.64| ) was evaluated using the lowest order profile for the domain 



wall, 

h(x x ) = V - [tanh(mx ± /2) + 1] , (3.85) 

with v being the true (not just leading order) scalar vacuum expectation value in the Higgs 
phase. Corrections to the profile of the domain wall only affect a at next-to-next to lead- 
ing order in the small 77 expansion ( |3.84| ). (This is tantamount to using the unperturbed 



wavefunction to compute the first order correction to the energy of a state in quantum 
mechanics.) 

^From Fig. p] it is clear that the difference between the second order derivative expansion 
result (|3.80|) and the full one- loop result is quite small, at least for X/g 2 < 1. Numerically, 
this difference appears to scale as 77 s / 2 . In contrast, the "classical" result ( |3.70| ) which omits 
all one-loop derivative corrections, and the leading two terms of the asymptotic expansion 
( j3.84| ), both depart significantly from the full one-loop result for quite small values of r\. 



IV. CONCLUSION 

The method presented in this paper provides a straightforward technique for evaluating 
the one-loop energy density of any configuration which depends non-trivially on only a single 
coordinate. It allows ultraviolet divergent contributions to be clearly isolated and extracted 
using any convenient regulator (which does not constrain momenta in the one non-trivial 
direction). As shown by our example of scalar QED, the method remains applicable in 
gauge theories, where multiple fluctuating fields are unavoidably coupled. Generalizing the 
treatment to, for example, full electroweak theory is straightforward. 

Note added: After the completion of this work, we became aware of Refs. P^T9|, in 
which methods quite similar to those we describe are applied to the closely related problem 
of fluctuation corrections in bubble nucleation. 
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APPENDIX A: FUNCTIONAL DETERMINANTS FROM THE SOLUTIONS OF 
ASSOCIATED DIFFERENTIAL EQUATIONS 

The ratio of one-dimensional determinants ( |2.32| ) follows from the result that, for any 
bounded functions U^°'(x±) and U^'(x±), 



det ± (-dl + U( \x ± )-X) _^(X;L) 



(Al) 



det ± (-9i + f/( 1 )(x ± )-A) <f(X; L) ' 
where ip an d f are the solutions to the associated ordinary differential equations, 

~-d\ + f/ (0) (x x ) - X] i/j(X; x ± )=0, (A2) 
~-dl + U (1 \x ± ) - X\ V {X;x ± ) =0, (A3) 
with boundary conditions specifying vanishing initial values 

</>(A; —L) = 0; <p(X; —L) = 0, (A4) 

and unit initial slopes, 

V(A; -£)' = !; (p(X; —L)' — 1. (A5) 



Here primes denote d/dx±, and the functional determinants in (|A1|) are defined on the space 
H.l of functions vanishing at x± = ±L. 

To prove the result (|A~1|), one need only notice that both sides of this equation are 
meromorphic functions of A with the same poles and zeroes, and that both sides approach 
one as |A| — > oo in any direction away from the positive real axis. Consequently they must 
be identical. This proof is due to Coleman (see Ref. Hl3| ). It may be easily generalized to 
the case of higher- dimensional operators. Let = ||W^° ||, a = 0, 1, be two self-adjoint 



second order differential operators acting on i= i.. n 7^L of the form = —dj_ + (x±). 
Then by the same argument sketched above, 

deta^W-AI) _ E n .., n e ll ... in 4 1 1 \^L)...4: ) (X;L) 



det ± (WW - XI) Eji .., neil .., ra ^)(A;L)...^(A ; L) 



(A6) 



where / is the unit matrix, and ^'s and <£>'s are the solutions of the associated system of 
ODEs, 

Y / (wS ) -X6 lJ )4 k \X-x ± )=0, (A7) 

j 

J2(W t f-X5 lJ )<pf\X;x ± )=0, (A8) 
j 

with the initial conditions 

# ) (A;-L) = 0; 4 k) (X; -L)' = 6 tk , (A9) 
^ (fc) (A; —L) = 0; #(A; —L)' = 8 ik . (A10) 
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The solutions {ip^(x±), ip( n \x±)} are linearly independent, and therefore the numera- 
tor in the right-hand side of ( |A6| ) vanishes if and only if there is a non-trivial solution of 
W^ijj(X; x±) = A'0(A;x_l) which vanishes at x± = ±L. In other words, the numerators of 
both the left and right hand sides of Eq. ( |A6| ) vanish whenever A is an eigenvalue of the 
linear operator W^ 1 ' (defined on the interval [—L,L] with Dirichlet boundary conditions). 
The same argument applies to the denominators on either side. One may also show that 
both sides approach unity as |A| — ► oo in any direction away from the positive real axis, and 
therefore conclude that both sides must be equal. 

The specific relation ( |3.44| ) immediately follows as the special case of n = 2, = A^, 
W ( ) = ^{o) ? and A = o. 
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